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Abstract 

We calculate properties of neutron stars such as mass and radius using a 

relativistic Dirac-Brueckner-Hartree-Fock approach. Modern meson-exchange 

potential models are used to evaluate the G-matrix for asymmetric nuclear 

matter. For pure neutron matter we find the maximum mass to be Af max ~ 

2AMq for a radius R ~ 12 km, whereas with a proton fraction of 30%, we find 

Af max ~ 2.1 Af© for a radius R ~ 10.5 km, close to the experimental values. 

The implications are discussed. 

PACS numbers: 97.60. Jd, 21.60.-n, 21.65. +f, 
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The properties of neutron stars depend on the equation of state (EOS) at densities up 
to an order of magnitude higher than those observed in ordinary nuclei. Data on the EOS 
can be obtained from many sources, such as studies of the monopole resonance in finite 
nuclei, high-energy nuclear collisions, supernovae and neutron stars. Supernova simulations 
seem to require an EOS which is too soft to support some observed masses of neutron stars, 
whereas analyses of high-energy nuclear collisions indicate a rather stiff EOS, predicting 
neutron star masses which are too large. Thus, no definite statements can be made about 
the EOS at high densities, except that it should probably be moderately stiff in order to 
support maximum neutron star masses in a range from approximately 1.4M Q to 1.9M Q [|TJ. 

However, although quantitative calculations of the EOS for dense nuclear matter are still 
beset with many problems, there have recently been important changes in the qualitative 
picture. Several mechanisms have been studied, mechanisms which result in an EOS which 
is soft enough to support neutron stars with maximum masses in the range of the observed 
data. Among such processes we find exotic states of nuclear matter, such as kaon [|2]|| 
or pion condensation ||. Another scenario which gives neutron star masses within the 
experimental values, has been presented by Pethick and co-workers These authors 

study properties of various phase transitions from spherical nuclei to uniform nuclear matter, 
and the coexistence of quark matter and nuclear matter over a finite fraction of the neutron 
star volume. 

The scope of this work is to derive the EOS for asymmetric nuclear matter, using rel- 
ativistic many-body theories within the framework of the Dirac-Brueckner-Hartree-Fock 
approach [|7||| and treating the Pauli operator which enters our formalism (see below) cor- 
rectly. To our knowledge, this has not been done before. Asymmetric nuclear matter is 
important in e.g. studies of neutron star cooling, as demonstrated recently by Lattimer et 
al. Q who showed that ordinary nuclear matter with a small asymmetry parameter can cool 
by the so-called direct URCA process even more rapidly than matter in an exotic state. In 
addition, at high densities, degrees of freedom represented by isobars and hyperons may also 
result in a rapid cooling of neutron stars [|10] . 



Before we present our results, we briefly sketch below our calculational procedure, for 



further details see Ref. [|TTJ. Important ingredients in our calculations are the nucleon- 
nucleon (NN) interaction which we take from meson-exchange models, and the renormalized 
NN potential in a medium like nuclear matter. The latter is accounted for by the reaction 
matrix G 

GM = V + VQ l—QGiu), (1) 

to — C^-Hov 

where u is the energy of the interacting nucleons, V is the free NN potential, Hq is the 
unperturbed energy of the intermediate scattering states, and Q is the Pauli operator which 
prevents scattering into occupied states. The G-matrix is used in our many-body scheme 
described below in order to obtain the energy per particle in asymmetric nuclear matter. 
To calculate the NN potential, we use the Bonn A potential defined in table A.l of Ref. 



12|| . This potential was shown by Brockmann and Machleidt |8j, to reproduce the nuclear 
matter binding energy and saturation density within the framework of the Dirac-Brueckner- 
Hartree-Fock approach. 

In order to test whether our EOS is appropriate for neutron stars, we evaluate the total 
mass and radius of such stars. Although there are recent relativistic calculations which 



address asymmetric nuclear matter, see e.g. the recent work of Huber et al. ||13|| , one of 
the differences between this work and Ref. [[HJ is an exact treatment of the Pauli operator 
Q for asymmetric nuclear matter. The conventional approach has been to extrapolate the 
results for symmetric nuclear matter and pure neutron matter ||14|| . Here we employ an exact 



expression for Q in nuclear matter as described in e.g. Ref. |15| and calculate the energy per 
particle for asymmetric nuclear matter using an extended Brueckner-Hartree-Fock (BHF) 
method, see Ref. [16j for a discussion of the non-relativistic BHF theory, namely the Dirac- 
Brueckner-Hartree-Fock (DBHF) method. The DBHF method is a variational procedure, 
where the single-particle (sp) energies are obtained through an iterative self-consistency 
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scheme. To describe the sp properties we depart from the Dirac equation for a free nucleonf], 



i.e., 



(i /) — m)ip(x) = 0, 



(2) 



where m is the free nucleon mass and ip(x) is the nucleon field operator, which is conven- 
tionally expanded in terms of plane wave states and the Dirac spinors u(p,s), and v(p,s), 
where p = (p°, p) is a four momentum and s is the spin projection. The relativistic energy 



is E(p) = \Jm? + |p| 2 . To account for medium modifications to the free Dirac equation, we 
introduce the notion of the self-energy E(p), which for nucleons can be written as 



E(p) = E s (p)- 7o E (p)+ 7 pS». 



(3) 



The momentum dependence of E° and Eg is rather weak |||17]. Moreover, E y << 1, and 
it is customary to rewrite the self-energy as 



E^E s - 7o E = f/ s + £/ v 



(4) 



where Us is an attractive scalar field and Uy is the time-like component of a repulsive vector 
field. The finite self-energy modifies the free Dirac spinors as 

/ \ 



u(p, s) 



' E{p) + m 
2m 



Xs 



a p 



(5) 



where we let the terms with tilde represent the medium modified quantities. Here we have 
defined [0,0 



m = m + Us 



(6) 



and 



1 Hereafter we set G = c = fr = 1, where G is the gravitational constant. 
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E a = E(p a ) = ^m* + pl (7) 
and the sp potential is given by 



~ 1 

22 -~-~-(ah\G(u = e a + e h )\ah) AS , (8) 

h<k F Eh^a 

where G is the relativistic G-matrix ||. We can also express the sp potential in terms of 
the constants Us and Uy, 

Tfl 

u a = -~rU s + Uy. (9) 
The sp energies e can then be written as 

e a = E a + U v , (10) 



where we have used the continuous choice for the single-particle energies [jl6f| . 

Eqs. are solved self-consistently starting with adequate values for the scalar and 

vector components Us and Uy. This iterative scheme is continued until these parameters 
show little variation. The calculations are carried out in the nuclear matter rest frame, 
avoiding thereby a cumbersome transformation between the two-nucleon center-of-mass sys- 
tem and the nuclear matter rest frame. The additional factors ih/E in the above equations 
arise due to the normalization of the neutron matter spinors w, i.e. w*w = 1 ||. With 
the obtained sp energies, we can calculate the relativistic energy per particle, see e.g. Refs. 



|8|JTT[| . All relevant equations with a given proton- neutron fraction are given in Ref . [|T^ . 

With these preliminaries, we present our results for the energy per particle with various 
proton fractions in Fig. |l} For relatively small proton fractions, the energy per particle 
exhibits much the same curvature as the curve for pure neutron matter at high densities, 
although the energy per particle is less repulsive at high densities. At lower densities, the 
situation is rather different. This is due to the contributions from various isospin T = 
partial waves, especially the contribution from the 3 Si- 3 Di channel, where the tensor force 
component of the nucleon-nucleon potential provides additional binding. From Fig. |T] we 
note that with a proton fraction of 15%, the energy per particle starts to become attractive 



at low densities (in the region 0.07 fm -3 to 0.3 fm -3 ). For larger proton fractions, additional 
attraction to the energy per particle is introduced. The next step in our calculations is thus 
to evaluate the EOS and the total mass and radius of a neutron star from the above energies 
per particle with the proton fractions shown in Fig. [I], in order to see how different proton 
fractions influence the mass and radius of neutron stars. Here we assume that the neutron 
stars we study exhibit an isotropic mass distribution. Hence, from the general theory of 
relativity, the structure of a neutron star is determined through the Tolman-Oppenheimer- 
Volkov equation, i.e. 

dP _ _{p(r) + P(r)}{M(r) + 4nr 3 P(r)} 

~~dr ~ r 2 _ 2rM(r) ' ^ ' 

where P(r) is the pressure and M(r) is the gravitional mass inside a radius r. To obtain 
observables like the mass and radius of a neutron star, we combine Eq. ([11]) with the equation 
of state (EOS), which is defined as P(n) = n 2 (de/dn), where e = S/A is the energy per 
particle and n is the particle density. Our EOS is valid in a limited density range from 0.1 
fm -3 to 0.8 fm -3 . It is therefore coupled to other equations of state at higher and lower 
densities as outlined in Ref. jll|| . Total masses and radii are calculated and parametrized 
as functions of the central density n c . These results are presented in Figs. and ^, for the 
total mass and radius, respectively. 

The neutron star equation of state should probably be only moderately stiff to support 
maximum neutron star masses of only 1.9M 0. From Figs. ^] and |] we see that our 
relativistic EOS for pure neutron matter seems to be too stiff, since it gives a predicted 
maximum mass of M max m 2.4M Q with a corresponding radius of R = 12 km. However, the 
EOS for neutron star matter could be softened considerably due to pion or kaon condensation 
because of the lower symmetry energy of nuclear matter, and maximum masses are then 
reduced correspondingly from the cases with no condensates. We see from Figs. ^] and |3| 
that for our EOS the calculated maximum mass can be reduced to M max w 2.OM with a 
corresponding radius of R ~ 10 km. Pions may be likely to condense in neutron star matter 
because neutrons at the top of the Fermi sea could decay to protons plus electrons f|, and 



kaon condensation is also believed to be a possible mechanism which could be energetically 
favorable in the interior of neutron stars P JI8| , |19|| . Both pion and kaon condensation would 
then increase the proton abundance in the matter, possibly up to more than 40% protons, 
i.e. close to symmetric nuclear matter, and produce a softer EOS and smaller maximum 
mass. Our results show this, if we assume that e.g. kaon condensation is a likely mechanism, 
an increased proton fraction results in smaller masses. Our EOS, even with a proton fraction 
close to symmetric nuclear matter, results in maximum masses which are slightly above the 
experimental values Q. However, our calculation of the EOS is to first order in the reaction 
matrix G, and we would therefore expect that higher-order many-body contributions to 
soften the EOS further. This was indeed shown in a preliminary study for symmetric nuclear 
matter by Jiang et al. [pOR . Although only a set of higher-order contributions was considered, 
the above authors obtain a softening of the relativistic EOS. The study of such effects will 
be studied by us in future works. 
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FIGURES 



FIG. 1. The energy per particle for asymmetric nuclear matter as function of the particle 
density for different proton fractions. 
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FIG. 2. M/Mq for various proton fractions as function of the central density n c . 
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FIG. 3. The total radius R for different proton fractions as function of the central density n c . 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9406028vl 
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